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*vj ' Abstract. We constructed Hirota-Kimura type discretization of the 

classical nonholonomic Suslov problem of motion of rigid body fixed 
at a point. We found a first integral proving integrability. Also, we 
have shown that discrete trajectories asymptotically tend to a line of 

(yr\ , discrete analogies of so-called steady-state rotations. The last property 

completely corresponds to well-known property of the continuous Suslov 
case. The explicite formulae for solutions are given. In n-dimensional 

1^"! , case we give discrete equations. 
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1. Introduction 



There are several methods for constructing discrete counterparts of infe- 
ct ' grable dynamical systems in classical mechanics. One is well-known Veselov- 
Moser discretization (see [12] )• By this method, based on discrete variational 
principle, many of discrete integrable systems are found (see [H])- In cases 
with discrete Lagrangian, corresponding discrete map is Poisson with respect 
to certain Poisson structure. Usually the discrete map is multi-valued. Re- 
cently, Hirota and Kimura constructed explicit integrable discretizations of 
the Euler and the Lagrange cases of motion of a heavy rigid body fixed at 
a point using Hirota's bilinear method [51 [9] . They found first integrals of 
motion and they solved equations in terms of elliptic functions. Suris and 
Petrera in [13] found bi-Hamiltonian structure for discrete Euler top. 

Our goal is to apply Hirota-Kimura method to the classical nonholonomic 
Suslov problem. Suslov in [15j considered motion of a rigid body fixed at a 
point with projection of angular velocity to an axis fixed in the body equal 
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to zero. Solutions in terms of trigonometric and exponential functions are 
given as well in [15]. There are various generalizations of classical Suslov 
problem: integrable potential perturbations ([10112]) and higher-dimensional 
generalizations ( [31 [8l [71 [16]). In ^ Fedorov and Zenkov presented certain 
discretization of Suslov problem based on the Veselov-Moser discretization 
and its extension to nonholonomic cases suggested in [H |TT]. Obtained 
map is multi-valued. Discrete trajectories asymptotically tend to discrete 
analogies of so-called steady-state rotations, the property characteristic for 
continuous case considered by Suslov. Moreover, Fedorov and Zenkov [4] 
gave the equations in n-dimensional case and proved similar asymptotic 
behavior as in three-dimensional case. 

In the present paper we use the Hirota-Kimura method and get explicit 
discretization of the reduced Suslov problem on the linear subspace of alge- 
bra so(3). The corresponding map is not multi- valued. We gave one first 
integral which appears to be enough for integration. Using a linear change 
of variables, we have found explicite solutions in terms of exponentional 
functions. (Notice, that in the case of discretization given in [1], the ex- 
plicite solutions are not known.) Presented discrete version of Suslov case 
has similar asymptotic behavior as the continuous one. 

The paper is organized as follows. In Section 2 basic facts about the 
classical Suslov problem are given. In Section 3 we present discrete equations 
of Hirota-Kimura type and we construct first integral of motion. Integration 
procedure for discrete equations is performed in Section 4. In Section 5 we 
give the discrete equations for n-dimensional Suslov case. 



2. A BRIEF ACCOUNT OF THE CLASSICAL SuSLOV PROBLEM 

The classical nonholonomic Suslov problem is defined in ^15j. Configura- 
tion space is Lie group 5*0(3). In a bases fixed in the body, the equations 
of the motion are: 



M = M X n -F Aa 

(1) 

^ ' (a, U) = 0. 



Here il is the angular velocity, and M = /i7 is the angular momentum, 
/ is the inertia operator, a is a unit vector fixed in the body and A is the 
Lagrange multiplier. In a bases chosen such that a = (0, 0, 1) and 
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(2) 



the equations ^ become: 

h^l = — -^13f^l^2 — -^23^2 
12^22 = hs^l + 123^1^2 

= -hsni - 123(^2 + {h - 12)^1^2 + A 

^3 = 0. 

The first two equations are closed in f^i and i^2- After solving them 
one finds the Lagrange multiplier A as a function of time from the third 
equation. Hence, for complete integrability by quadratures, only one first 
integral of motion is necessary. This is the energy integral as follows from ([2]) 
easily. Suslov in [15] gave solutions of the system in terms of trigonometric 
and exponential functions. He observed a remarkable fact (as the referee 
observed to be known before to Walker and Routh): the motion of the body 
asymptotically tends to a line of rotations with constant angular velocities 
which satisfy /is^Ji + /23^2 = 0. 

3. HiROTA-KlMURA TYPE DISCRETIZATION OF THE SuSLOV PROBLEM 

In the spirit of Hirota-Kimura, a discrete counterpart of the first two 
equations and nonholonomic constraint of ([T]) is: 






^(n^n^ + Q^ns) 



+ /23f^3!^3 



-(02^^3 + 02^3 



'13 



(O1O2 + O1O2) - /23f^2^2 



(3) 



l2(02-02)+/23(f^3-^3 



-(O1O3 + O1O3 



O3 



/13O3O3 - Y(f^lf^3 + ^lf^3) + ^(^lf^2 + f^lf^2) + /l3^lS^l 
^ -O3. 



Here Oj = Oj(t), Oj = Oj(t + e) and e is the time step. The limit when e 
goes to should reconstruct equations ([3]) of the continuous Suslov problem 
([2|). Thus one concludes that O3 should be equal to zero, and equations ([3|) 
become: 

^3, 



(4) 



/i(Oi -Oi) = e 

12(^2 -n2) = e 

03 = 



-(Oi02 + Oi02)-l23^2f^2 



/. 



23 , 



^-^{^1^2 + ^1^2) +113^1^1 



Since these equations are linear in Oj, the map defined by (jl]) is explicit and 
unique- valued: 
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giving 



(5) 



where 
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^ A^ 2/2 2/1 h ^' 

A ^ 2/1 2 2/2 /2 



As we have aheady mentioned, in continuous case equations ([2]) have the 
energy integral. But for considered discretization the energy is not an in- 
tegral anymore, but there exists a first integral as it can be seen from the 
following statement. 

Lemma 1. The function 

hnl + hnl 

^^' "4Iil2+e2(/i3fli+/23f^2)2 

is a first integral of equations ([5|) 

Proof. Proof follows by direct calculations. D 

In the limit when e goes to zero, integral ([6]) tends to the energy integral 
divided by constant. 

4. Integration 

In order to integrate discrete Suslov equations, we introduce new coordi- 
nates: 

(7) X = /isili + 723^2, y = 123^^1 - 713/2^2- 
The Jacobian of the change of coordinates ([7]) is: 

-iihh + ihh). 

Thus, it is equal to zero only in the case /13 = /23 = when ill = Qi and 
^2 = ^2; giving equilibrium position. 

In the new coordinates d?]) equations (JH) are: 

x-x= ^rr-j-ixy + xy) 

(8) ^^1^2 

y — y = —exx. 
The first integral ([6]) becomes 

(9) p^hl2X^ + ll_ 

^^ 4Iil2 + e2x2- 
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The curve F{x, y) = h can be parameterized by introducing: 



hhh 

^ = ^\ Tl 1 2^^°*^' 

(10) ^ ^ 



h 
y = 2Iil2\ —- 2:;rsin0. 

Putting (fTOl) into the second equation of ([8|) and denoting 



hh 



we get: 



(11) nv^yTT-s7^:^TT= ,%"^, - 

Let us suppose the form of a solution of ([TT]) : M(n) = sh(/i:i(to + f^e) + ^2)1 
where ki and A;2 are constants. By plugging the form into (jlip . one gets: 

(12) sh(-feie) = jj^^j^,- 

So, we have 

Proposition 1. Lei constant ki satisfy (I12p anrf A;2 6e arbitrary constant. 
Then the function u{n) = sh(/ci(to + ?^e) + ^^2) gives solutions of equation 

m- 

As a consequence we have the following statement: 

Theorem 1. Let constant ki satisfy (jl2[) and k2 be arbitrary constant. The 
functions: 

2Vh / /13/2 

(Ifg/i + If^h) ch(/ci(io + ne) + fez) V Vhh^^W 

2^/h / hzh 

(Ifg/i + /fg/a) ch(A:i(fo + ne) + ^2) V ^hh - he' 
give the solutions of equations ([5]) . 

Proof. From Proposition [1] and (jlOp one gets: 



^i("') = /r2 r _L r2 r \ i.^> /^ — , ^ \ , i \ \ It t , 7 + -^23 sh(A;i(to + ne) + k2, 
IS-*! + -'13 " " 

^2('^) = ij2 J , t2 T \^^n. u — , „,A , ,„ A ( ^T=^=^=7 ~ ^13 sh(A:i(to + ne) + ^2) 



, . „ / /l/l^2 1 / \ o /7n-l~^^(^l(*0 + '^'^) + ^2) 

^^ V /i/2-/ie2ch(A;i(to + ne) + A;2)^^ ^ ^ ' ch(A;i(to + ne) + ^2) 
Proof follows from ([7]) . D 

Proposition 2. The discrete trajectories of motion of the body asymptot- 
ically tend to a line of discrete analogies of so-called steady-state rotations 
that satisfy: 

hsfti + 123^2 = 0. 
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Proof. In the limit n — > ±00, we have that x goes to zero giving that 
Iiaili + /23^2 goes to zero. D 



The last statement is illustrated by the following pictures. In all four 
cases, the line /13O1 + 123^2 = is clearly indicated. Following parameters 
are chosen for picture 1: e = 0.2, Ii = 4, I2 = l,/i3 = — 0.5,l23 = —0.3, 
for picture 2: e = 0.2, h = 4, /2 = 3, In = -0.4, I23 = -0.2, for picture 3: 
e = 0.02, Ii = 4,/2 = 2,Ii3 = 0,/23 = -0.2 and for picture 4: e = l,Ii = 
3,l2 = 3,/i3 = -0.2,/23 = -0.2. 
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Remark 1. Equation (|lip splits on two equations: 



u 



zy V? + 1 + uy (? + 1, 



-cyv?' + 1 — uy (? + 1 



where c = j H ^^'i ■ The solutions given in Proposition [T] correspond to the 
first equation. One easily concludes that for the second equation the limit 
when e goes to zero is not defined well. 

Remark 2. One can use the change of coordinates x = Irs^i + 123^2, y = 
l23li^i — Ii3l2^2 also in the continuous case. Then equations ^ become: 



xy 

hi 2 



y 



-x\ 



which are simpler then original ones. One can easily see that Hirota-Kimura 
type discretizations of last equations are equations ([7]) , as linear change of 
variables commutes with Hirota-Kimura type discretization. 
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5. Higher-dimensional case 

As it has already been mentioned, the higher-dimensional generalization 
of Suslov case was suggested by Kozlov and Fedorov (seeOHj). The config- 
uration space is Lie group SO{n) and the nonholonomic constraints are: 

i^ij = 0, 1 < i,j < n — 1. 

The equations of motion are: 

(13) M = [M, f)] + A 

where M = Ifl + ClI, and 



hi ... hn 
I22 ... hn 



Lin J-2n 
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A 




— A12 

-Al,n-1 




A12 



-A2,n-1 





Al,n-1 
A2,n-1 






From (|13p closed systems of equations in VLm-, 1 < i < n — 1 follows: 

(14) 

y-'-ii ~r -InnJ^Hn ^ -'inv^ 'in ' ••• ~r ^^n—l,n) ~T~\'^ln^^ln + ... + -'n— l,n" ^n— l,raj" H' 

Similarly as in three-dimensional case, we give the Hirota-Kimura dis- 
cretization in n dimensions by the following system of equations: 



\-i a ~T~ ^ nn ) \^ Hn ^Hr, 



(15) 



— ^^in\^^ln^^ln i ■■■ + ^^n~l,n^^n— l,n 
+ ^{hn^ln + ... + h-l,n^n-l,n)—:r- 

+ ^{hnP'ln + ... + In-l,n^n-l,n)—;z-- 



As in three-dimensional case, map defined by ()15p is explicit and unique- 
valued: 



(16) 



where 



A 



^m 


= A-' 


' f^m " 


.f^n-l,n. 




"n— l,n 



e{l2n^2n + ■■■ + !■, 



n — l,n^^n — l,n) 



2(/ll+/nn) 
e{2l2n^ln—h,n^2n) 



'2{l22+lnn) 



^\-tn—l,n — l~r-L'nn) 



2(/ll +/nn) 

e(2J2n^n-l,n— -fn-l.n^2n) 

2 (/22+-f nn) 



1 £(-fln^ln + ...+Jn — 2,ra^ra — 2,n) 

^\-tn — l,n—l~r-Lnn) 
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Let US mention that equations ()15p have a particular solution -j^ = ... = 
j"~^'" = c = const, which corresponds to motion with constant angular 
velocity. 
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